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Abstract. We consider the inverse problem of the determining the potential in the 
dynamical Schrodinger equation on the interval by the measurement on the whole 
boundary. Provided that source is generic using the Boundary Control method we 
recover the spectrum of the problem from the observation at either left or right end 
points. Using the specificity of the one-dimensional situation we recover the spectral 
function, reducing the problem to the classical one which could be treated by known 
methods. We adapt the algorithm to the situation when only the finite number of 
eigenvalues are known and provide the result on the convergence of the method. 

1. Introduction 

We consider the problem of determining the potential for a one dimensional Schrodinger 
equation from two boundary measurements. More precisely, given q £ L^(0,1) (without 
loss of generality, we assume that q is real- valued) and a G Hq{0,1), we consider the 
following initial boundary value problem: 

iut{x, t) — Uxxix, t) + q{x)u{x, t) = t > 0, 0<x<l 

u{0,t) =u{l,t) = t>0, (1.1) 

u{x,0) = a{x) < x < 1. 

Assuming that the initial data a is unknown, the inverse problem we are interested in is 
to determine the potential q from the trace of the derivative of the solution u to (II. ip on 
the boundary: 

{ro{t),ri{t)} :={n,(0,t),n,(l,t)}, Vt E (0,r), 

where T > is a fixed constant. Once the potential has been determined, we can use the 
method of iterative observers recently proposed in [16j to recover the initial data. 

The multidimensional inverse problems for the Schrodinger equation by one measure- 
ment were considered in [41, [5], (Si I14j . Using techniques based on Carleman estimates, 
the authors established global uniqueness and stability results for different geometrical 
conditions on the domain. We also point out the approach developed by Boumenir and 
Tuan for the heat equation in [9l I10|, 111] . Using observation at the one end of the interval 
for t £ (0,oo), the authors were able to recover the spectrum of the string. Next, choos- 
ing another boundary condition, they solved the inverse problem from the two recovered 
spectra by the classical Levitan-Gasymov method jl2^ I13j . 

In this paper we offer a different approach. More precisely, let us introduce the un- 
bounded operator A on L^(0, 1) defined by 

Acp := -(I)" + qcj), y (p eV{A) := H'^{0,l)nH^{0,l). (1.2) 



Using the Boundary Control Method (denoted BCM in the rest of the paper, see [U [6]), 
we first show that the eigenvalues of A can be recovered from the data tq (or ri). To 
achieve this, we derive a generalized eigenvalue problem involving an integral operator 
(see (j2.5p ). whose resolution leads to the recovery of the eigenvalues of A. Then, using 
the peculiarity of the one dimensional case, we recover the spectral function associated 
with A, reducing the initial inverse problem to the more "classical" one of recovering an 
unknown potential from spectral data. We can then use either Gelfand-Levitan, Krein or 
Boundary Control Methods (see \2\). 

Remark 1.1. The method of solving the inverse problem presented in this paper could be 
applied to the case of wave and parabolic equations on the interval as well. 

The paper is organized as follows : in Section O we detail the different steps of our 
algorithm. In particular, we describe how to recover the spectral data of A from the 
boundary data and point out two methods that can be used to recover the potential. In 
Section [31 we show how to adapt our algorithm to the more realistic situation where not 
all but only a finite number of eigenvalues of A are known. 



2. Inverse problem, application of the BCM 

2.1. From boundary data to spectral data. It is well known that the selfadjoint 
operator A defined by (II. 2p admits a family of eigenfunctions {4'k}^=i, associated to a 
sequence of eigenvalues — >• +oo, forming a orthonormal basis of L^{0,1). Using the 
method of moments, we can represent the solution of (11. ip in the form 

oo 

uix,t) = Y,ake'^'''Mx), (2.1) 
1 

where are (unknown) Fourier coefficients: 

ak = {a,(l)k)L^o,i;dx)- (2-2) 
Definition 2.1. We call the initial data a G i^o(0, 1) generic if Uk for all k > 1. 

For the boundary data tq, ri, we have the representation 

{oo oo "\ 

ake^'-'4>',{0), ^ke^'^'^'kil) • (2.3) 
k=l k=l ) 

In particular, one can note that ro,ri G L^(0,T). Throughout the paper, we always 
assume that the source a is generic. The reason for this requirement can be observed 
from (12. 3p : if we assume that a„ = for some n, then from the representation (12. 3p it is 
easy to see that the response function does not contain any information about the triplet 
{An, 0^(0), (/>^(1)}, and, consequently the potential which would be possible to recover 
from {ro(t), ri(t)} would not be unique. 

Using the method from [l] we can recover the spectrum Afc and products afc(/)'^(0) and 
Ofc</'fc(l) by the following procedure. We construct the operator : L^(0,T) — > L^(0,T) 
by the rule: 

iC^)fit)= r roi2T-t-T)fiT)dT, OsSt^T. (2.4) 
Jo 
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and consider the following generalized eigenvalue problem : Find /) G C x L^(0,T), 
/ / 0, such that 

r ro{2T-t-T)f{T) = r ro{2T-t-T)f{T)dT, O^t^T. (2.5) 
Jo Jo 

Then, one can prove that the above problem admits a countable set of solutions {fin, fn), 
n ^ 1. Moreover, the eigenvalues /i„ are precisely given by the values iXn, where are 
the eigenvalues of A and the family {/n(i)}^i is biorthogonal to {e*'^'=*^-^~*)}^^. 

Next, we consider the equation of the form (|2.5p with ro(t) replaced by its complex 
conjugate ro(t). This equation yields the sequence {—i^k, 9k(t)}'^^^, {—iXn = /«„)• Let us 
normalize functions fk, Qk by the rule: 

Snk = {C^fn,9k) = / ro{2T-t-T)UT)gk{t)dTdt (2.6) 
Jo Jo 

and introduce constants 7^, (3k defined by: 

7fc= r ro{T-T)fk{T)dT, (2.7) 
Jo 

(3k = [ ro{T-T)'^dT. (2.8) 



Jo 

Then the product akcp'^iO) could be evaluated by 

ak(()'k{0) = IkPk. (2.9) 

Similarly we can introduce the integral operator Cj associated with the response at the 
right endpoint ri(i), and repeat the procedure described above to find quantities afc(/>'^(l). 

Summing up, using the method from [Ij we are able to recover the eigenvalues \k of 
A and the products (/)'^(0)afc and 0'^(l)afc. As a result, we can say that we recovered the 
spectral data consisting of 

Precisely this data was used in |15j . where the authors shown the uniqueness for the 
inverse spectral problem and provided the method of recovering the potential. Instead of 
doing this, we recover the spectral function associated to A and thus reduce the inverse 
source problem to the classical one of determining an unknown potential from spectral 
data. 

Given A G C, we introduce the solution y(-. A) of the following Cauchy problem on (0, 1): 
-y"{x,\)+q{x)y{x,\) = \y{x,\), < x < 1, (2.11) 
y(0,A)=0, y'(0,A) = l. (2.12) 

Then the eigenvalues of the Dirichlet problem of A are exactly the zeroes of the function 
7/(1, A), while a family of normalized corresponding eigenfunctions is given by (/>fc(x) = 
y{x,\k) 



|y(-,A 



k) 



-. Thus we can rewrite the second components in T> in the following way: 

^^_2/(Mfc)_, 

Let us denote by dot the derivative with respect to A. We use the following fact (see [15] ). 
which is also valid for q G L^{0, 1): 
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Lemma 2.2. If Xn is an eigenvalue of A, then 

\\y{;Xn)\\l2=y'{l,Xn)y{l,Xn) (2.14) 

The lemma below which can be found in [15j for the case q G -^^(0, 1) holds true for 
q £ L^{0,1) as well. The meaning of it is that the function y(l, A) is completely determined 
by its zeroes, which are precisely the eigenvalues of A. 

Lemma 2.3. For q £ L^{0, 1) the following representations hold 

Afe — A 



2/(i,A)=n 



y{l,Xn) = ^ TT 



Xk — Xn 



Lemma 12.21 and 12.31 imply that the data D we have recovered (see (|2.10p ) allow us to 
evaluate the norm 

\\yi;Xn)\\l2=AnBn=:al (2.15) 



where we denoted 



Bn:=-^ n (2-16) 

k^l,k^n 

2.2. Reconstructing the potential from spectral data. The set of pairs 
{'^fci l|y(") ■^A-OllialfcLi is a "classical" spectral data. The potential can thus be recovered 
by Gelfand-Levitan, Krein method or the BCM (see [2]). Below we outline two possible 
methods of recovering the potential. 

We introduce the spectral function associated with A: 

E ^ A^o 

P(A) 



E 7^ A>0, (2.17) 



0<Afc^A 



which is a monotone increasing function having jumps at the points of the Dirichlet spectra. 
The regularized spectral function is introduced by 

0<A^ ^A 

where po is the spectral function associated with the operator A when q = 0. In the 
definition above eigenvalues and norming coefficients are A^ = vr^/c^^ (q^^)^ = 2T^k^ 

Sill V Ax 

the solution to dZTTTl - dm]) for g = is ^0(2;, A) = 

V A 

Let us fix r G (0, 1] and introduce the kernel c^{t, s) by the rule (see also [2]): 

/"°° sin vA(r — t) sin -v/A(t — s) , ,1 , 

c'{t,s)= — — -da{X), s,t£[0,T], (2.19) 

J-00 X 

Then so-called connecting operator (see [6l [2]) : L'^{0,t) i— )• L'^{0,t) is defined by the 
formula 

(CV)(t) = {! + nf{t) = fit) + r c^{t, s)fis) ds,0<t<T, (2.20) 

Jo 



Using the BCM leads to the following: for fixed r G (0, 1) one solves the equation 

iC-r){t) = T-t, 0<t<T, (2.21) 

Setting 

^(r) := r(+0), (2.22) 
and then varying r G (0, 1), the potential at the point r is recovered by 

We can also make use of the Gelfand-Levitan theory. According to this approach, for 
r = 1, the kernel satisfies the following integral equation with unknown V: 

Viy, t) + c^y, t) + r c^t, s)Viy, s)ds = 0, < y < t < 1. (2.24) 
Solving the equation (|2.24p for all y G (0, 1) we can recover the potential using 

q{y) = 2^V{T-y,T-y). (2.25) 
ay 

Once the potential has been found, we can recover the eigenfunctions (j)k, the traces 
4>'f,{0) and using (|2.9p . the Fourier coefficients a^, k = 1, . . . , oo. Thus, the initial state can 
be recovered via its Fourier series. We can also use the method of observers (see [16j). 

2.3. The algorithm. 

1) Take ro(t) := Ux{0,t) and solve the generalized spectral problem (j2.5|) . Denote 
the solution by {^ni/n (Ol^i note the connection with the spectra of A: 

2) Take the function ro(t) and repeat step one. This yields the sequence {ji^, 5n(i)}$S=i- 

3) Define the operator Cq by (12. 4p and normalize /„ , gn according to equation (12. 6|) : 

{Cjfn,gn) = l- 

4) Find the quantities a„(/)'„(0) by ([221), ([23]) and (l2:9]l . 

5) Repeat steps l)-4) for the function ri{t) := Ux{l,t) to find afci?!)'^(l). 

6) Define the spectral data T> by (|2.1Up and find the norming coefficients by using 
(f2T5]l and (f233|) . (f236|) . 

7) Introduce the spectral function p(A), the regularized spectral function cr(A) and 
the kernel respectively defined by (f2T7p . (f2T8]l and (f2T9]l . 

8) Solve the inverse problem by either BCM using ([MI]), (12221), or Gelfand- 
Levitan method using equations (|2.24p . (j2.25p . 

9) Use the method of iterative observers described in |16j or Fourier series to recover 
the initial data. 

Our approach yields the following uniqueness result for the inverse problem for 1-d 
Schrodinger equation: 

Theorem 2.4. Let the source a G -ffo(0, 1) be generic and T be an arbitrary positive 
number. Then the potential q G L^{0,1) and the initial data are uniquely determined by 
the observation {ux{0,t),Ux{l,t)} fort G (0,T). 

The method could be applied to the inverse problem for the wave and parabolic equa- 
tions with the potential on the interval as well. The details of the recovering the spectrum 
Afc and the quantities afc0'^(O), ak(p'j^{l) could be found in The following important re- 
mark, connected with the types of the controllability of the corresponding systems, should 
be taken into the consideration: 
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Remark 2.5. The time T of the observation could be arbitrary small for the case of 
Schrddinger and parabolic equations and is doubled the length of an interval (T = 2 in our 
case) for the wave equation with the potential. 



For the details see |T]. 



3. Stability of the scheme : the case of truncated spectral data 

In view of applications, in this section we consider the case where only a finite number 
of eigenvalues of A are available. More precisely, let us assume that we recovered the 
exact values of the first N eigenvalues A„ and traces An, for n = 1, . . . , N . Then we can 
introduce the approximate normalizing coefficients an,N by the rule 

an,N = AnBn,N, where Bn,N ■■= j-^ TT ^2 ^" (^•^) 



Then we can estimate 



1- n 

k^N+l 



(3.2) 



Since the infinite product ()2.16p converges, the right hand side of the above inequality 
is small enough as is getting bigger, provided n is fixed. But the following remark 
should be taken into the account. Let us remind the following asymptotic formulas for the 
eigenvalues and norming coefficients: 

\k = 'K'^k^ + c + o(^] , k^oo, (3.3) 



«l = ^+o(l), k^oo, (3.4) 

where c = q{s) ds. Then the product in the right hand side of (|3.2p can be rewritten 
Jo 

^ / n:^^n(X) / 1 \ \ 

(3.5) 





as 

n2+o(4) 

So we see that if n is close to N, then the terms j^^" ' are close to one, and consequently 

the first factors in (|3.5|) are small. This implies that the product in (|3.2p is not close to 
one. This simple observation yields that we can guarantee the good estimate in (13. 2p only 
if N is much greater than n, so that the product in (|3.5p is close to one. 



We need to find out how many eigenvalues (we call their number by N) we need to 
recover in order to n less than approximate normalizing coefficients be recovered with a 
good accuracy. In other words, we need to find out the relationship between and n such 
that the product in p. 20 or what is equivalent ()3.5p are close to one. Let us notice that if 
< X < ^ < 1 then | In (1 — x)| < Using this observation, we can estimate for 

6 



in n 1 



n 



^ |ln(l-g)| n^ _ 2I 



ln(l 



In ( 1 
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n 



^ 1 



fc2 



(3.6) 



fc^TV+l \ k=l 

If we choose N and n such that right hand side of (|3.6p is less than some e > 0, then 

" < 1. (3.7) 



ri' 
fc2 



n > 

Consequently, for such and n we have (see p.2p ): 

\an - an,N\ ^ |A„||5„,Ar| |l - e"^| . (3.8) 

Notice that since An = 1 + o(l) and a„ = ^ ^ + o(l) as n — oo, we have that |B„^7v| 
is bounded by some positive C < 4. So, taking e small enough, and picking and n such 
that the right hand side of p.6p be less that e, we arrive at ()3.8p . Using formula 



A;=l 



from [TTj we summarize all our observations in the lemma: 
Lemma 3.1. Let < e < 1. If n and N satisfy 

n2 



n 

— < 



e 



N ^ 2 In 2' 

then there exists an absolute constant C > such that 

\ak - ak,N\ ^ C |l - e""^! ,VA; = 1, 



(3.9) 



(3.10) 



Let us recall the following important representations (see |2j) for the response function 
and for the kernel c"^ (in our case r = 1): 

Lemma 3.2. Assume that q G L^(0, 1). Then the following representation for the response 
function r, 



r{t) 



sm 



V\t 



da{X), 



(3.11) 



holds for almost all t € (0, 2r). 

The kernel c'^{s,t) admits the following representation: 

sin \/A(r — t) sin ^/X{t — s 



c'{s,t) 



A 



da{X), s,t G [0,r] 



(3.12) 



where the integral in the right-hand side of Ii2.19\) converges uniformly on [0, r] x [0,r]. 



and p{t) is defined by 



c^(t, s) = p{2t -t- s)-p{t- s) 
1 fl*! 



p{t) 
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r (s) ds. 



(3.13) 
(3.14) 



The following useful formula follows directly from representations (j3.1ip . ()2.19p : 

1 /-SCr-i) 

c"(t,t)=2y^ r-(r)dr (3.15) 

If the exact values of the first n eigenvalues and normalizing factors a|, k = 1, . . . , n, 
were known, then one could construct the "restricted" response functions and kernels 
defined by 



k=l 



sin ^/X^t sign \l ^k^ 1 



6^0 {air 

K 



k=l 



sin ^/Xk{T — t) sin ^/Xk{T — s) sign Afc 



Afc 



at 



sin yAO(r - t) sin (T - s) ^ 



then every r„ G C°°(0, 1) and Lemma [32] yields 

r„(t) — ^ r(t), for almost ah t G (0, 2), 
c^(t, s) — )• c^(t, s), uniformly on (0, r)^, 
c^{t,t) — )• c^(t,t), uniformly on the diagonal. 



,0\2 ■ 



(3.16) 
(3.17) 
(3.18) 



As we only have access to approximate values of the normalizing factors in practice, we 
can only compute the approximate restricted kernel by 

~r I. ^ _ sin VA^(r - i) sin ^/A^(r - s) sign Afc 



k=\ 



Then we can estimate the difference 



Afc 



fe,JV 



sin ^l[T - t) sin ^A^ (T - s) ^ 



,0\2 ■ 



n 152 



fc,Af 



(3.19) 



(3.20) 



Using (j3.10p and asymptotical expansions for the eigenvalues and norming coefficients 
(13. 3p (13. 4p . we deduce from (I3.20p that (below, C denotes an absolute constant that might 
change from line to line): 



Afea^./vai 



■■'--^k,N^k 

Having remembered estimate ()3.10p . we finally get 



cYk\ 

k=l 



ak,N - ak\. 



,r ^ II ^ n(n + l) , I 



(3.21) 



Let us fix some 6 > and choose n G N such that ||c^ — c"^ 



difference 



I ^^"^7~ 'T' 1 1 ^ 1 1 "7" 1 1 1 1 'T' "7" I 

ICn.Af ~ C ||oo ^ l|Cn,A'' ~ l|oo + || C„ ~ C | 



^ — and consider the 



(3.22) 



Using (|3.2ip we obtain the existence of a constant C* > such that 

Wn,N - c'lloo ^ C*n' |l - e-^l + -. (3.23) 

Then by choosing e in ()3.23p we achieve 

\K,N-c^\\oo^5. (3.24) 

We summarize the above observations in the following Proposition, which details the main 
of our approximation procedure. 

Proposition 3.3. Let 6 > be fixed. Let n be chosen that (this is possible thanks to 

mi}) 
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Next, take e > such that 
Finally, choose N such that 



I t-n 1 1 oo ^ 2 ■ 



N 2 In 2 

Then, there exists an absolute constant C > such that the following estimate holds true: 

11^, Af ~ c'^lloo ^ C6, 

where the approximate restricted kernel c^ ^ and the approximate normalizing coefficients 
o.k,N, k = 1, . . . ,n are respectively defined by i3.19\} and ^3. 1]) . 

In particular, ^ converges uniformly to c^ when n tends to infinity and N is chosen 
as above. 

Along with the equation ()2.24p we consider the equation for the approximate restricted 
kernel c^^^: 

K,Af(y,t) + c^,Ar(y,t) + / c^^Nit,s)Vn,N{y,s)ds = 0, 0<y<t<T. (3.25) 

Picking 5 > we can use Proposition 13.31 to find n and such that ||cj^^ — £"^1100 ^ S. 
Prom now on, we always assume that n and N are chosen according to Proposition 13.31 
Note that in particular, we have N — )• +oo as n — )• +oo. 

We introduce the operator C"'^ which are given by an integral part of (j2.20p with the 
kernel c"^ substituted by ^. The closeness of c"^ and ^ implies the operator / + 
to be invertible along with I + Cr. The latter in turn implies the existence of the potential 
Qn N that produces the response function 



''n,Nit) = 



sm 



(3.26) 



k=l 

and the unique solvability of (f3:25|) (see [2], [7], [3]). 

We denote M := ||(/ + C^)-^||. The invertibility of I + and / + C"'^ implies the 
norms of the solutions \\V{y, and \\Vn^Niy, to be bounded. 
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Let us write down the difference ()2.24p and ()3.25p in tlie form 

V{y,t)-Vn,Niy^i)]+ [ c^(t,s) k(y,s)-K,7v(y,s)l ds (3.27) 



= c^,Af(y> t) - c^(2/> t)+ s) - Ar(t, s)] K,7v(y, s)ds, < y < t < T. 

The above equahty and the invertibihty of / + impUes the estimate: 



\\V{y, •) - K,7v(y, -^L^iy^r) ^M[l+ max \\Vn,Niy, OIIl-J I|c" " S^.Tvlloo- (3.28) 
To estimate the L°° norm, we write down (j3.27p in the form 

V{y, t) - VnMv^ t)] = r c^{t, s) \Vn,N{y, s) - Viy, s)] ds (3.29) 

+c^,N{y,t) - c'{y,t) + [c'{t,s) -c^^^{t,s)]Vn,Niy,s)ds, < y < t < T. 

Jy 

The latter leads to the following inequality: 

rT 

ds (3.30) 



Il^(y,*) - K,Af(y,t)||oo ^ llc^lloo max / Vn,N{y,s) -V{y,s) 



+ ||c^^ - c^lloo + llc"^ - c^atIIoo max / \Vn,N{y,s)\ds. 
Using ([3:28]) and (ICTjl we finally get 

llV'-K.Jvlloo ^ ||c^-2;jvlloo(M||c^||oo + l) (l + max \\Vn,N{y, 

The last equality in particular implies the uniform convergence of Vn,Niy,y) to V{y,y). 
Having remembered (|2.25p . we conclude that the potentials qn,N and q corresponding 
respectively to at and V satisfy 







/ q{s) ds, uniformly in t. 
Jo 



ln,N{s) ds =^n 

The latter in turn, implies that 

qn,N — > q, inif~^(o, 1). 

In fact, (I3.3ip yields more than that: we have 



t+e 



qn,Nis) ds 



1 r*+^ 

— / q{s) ds, uniformly in t, e. 
£ Jt 



(3.31) 
(3.32) 
(3.33) 



We remark that (I3.33P is still not enough to guarantee the convergence qn^N to q almost 
everywhere on (0, 1). 



On the other hand let us restrict ()2.24p to the diagonal y = t: 

rT 

V{y,y) + c^{y,y)+ c^{y,s)V{y,s)ds = 0, < y < 1, 

Jy 



(3.34) 



and recalling (j2.25p . ()3.15p . we see that the best possible result one can expect is a con- 
vergence qn^N to q almost everywhere on (0, 1). 
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Remark 3.4. The stability of the scheme crucially depends on the type of the convergence 
r„ — )• r. For now we know (see [2\) that the convergence is pointwise almost everywhere 
on the interval. The significant progress in the proving of the stability could be achieved 
by the improving of this result. 
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